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The NA48/2 collaboration has shown clear experimental evidence for a cusp in the data for K TTTr^Tr". This 
effect can be used to extract information on the nn scattering lengths. We address this issue using a two-loop 
dispersive construction of tttt — > tttt and K nmv amplitudes in the presence of isospin breaking. 



1. Introduction 

The observation of a cusp anomaly in the 
invariant mass distribution in the data on the 
K"^ TT^TrOTr" decay collected by the NA48/2 
collaboration [1] has triggered some theoretical 
activity [2]- [5]. The basic explanation of the ap- 
pearance of a unitarity cusp is very simple [2] : the 
amplitude for A'+ ^ tt+tt^tt" has two basic con- 
tributions, one of which corresponds to the tt+tt" 
intermediate state rcscattering to tt^tt*^. This in- 
termediate state in the s-channel generates (at 
the one-loop level) a square root singularity [6], 
and the corresponding amplitude behaves at this 
level for s ~ 4m^ (to+ is the mass of 7r+) as 

-4(.,.«(.i)+,sK,{i^-»/^t;.v w 

where ap = \/l ~ Arn^p/s, ap = -^/imlJ/s — 1- 
The functions i?(cr+) and S{a'^) can be expressed 
as convergent series in s — 4m^ in the physical 
region of the A'+ 7r+7r°7r° decay. This singu- 
larity appears at Arrv^, which is above the physi- 

*Presented by M. Z. and K. K. at Flavianct Kaon Work- 
shop 08 and QCD 08. This work was supported in part 
by the Center for Particle Physics (project no. LC 527), 
GACR (grant no. 202/07/P249) and by the EU Contract 
No. MRTN-CT-2006-035482, "FLAViAnef. 
In memoriam of our friend and colleague Jan Stern, 
t Unite Mixte de Recherche (UMR 6207) du CNRS et dcs 
Universites Aix-Marseille 1 et 2 et Sud Toulon-Var, La- 
boratoire affilic a la FRUMAM (FR 2291) 



cal threshold irriQ {mo is the mass of 7r°), and the 
cusp is a result of the interference of the part con- 
taining the singularity and the rest without it. It 
is clear that the cusp appears only in the isospin 
breaking case and its strength is sensitive to the 
TTTT scattering amplitude at the threshold. This 
is the reason why the investigation of the cusp 
effect can in principle serve as an independent 
method for the experimental determination of the 
TTTT scattering lengths, provided a model indepen- 
dent description of the corresponding amplitude 
can be given. 

In [3] Cabibbo and Isidori have proposed to use 
the assumed simple analytical properties of the 
amplitude (tttt amplitudes are considered in the 
form ([T]) with i?(o'^) and S{(j'^) polynomial) and 
the unitarity of the scattering matrix in order to 
express the amplitude near the threshold as an 
expansion in scattering lengths . This idea was 
further investigated in [?]. Another possibility 
how to address this issue is the framework of non- 
relativistic effective field theory as developed in 
[S]. This is done as a combined expansion again 
in Ui, and also in the pion momenta. Both of 
these different approaches compute contributions 
to the K'^ decay amplitude up to order 0{a1). 
Each of these approaches has its own limitations 
that result from the assumptions made. More- 
over, they give no connection between cusp effects 
and the traditional PDG parameterization of the 
K decay amplitudes. We consider an alternative 
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approach, which rests on general properties, uni- 
tarity, analyticity, crossing symmetry, relativistic 
invariance, and chiral power counting for partial 
wave amplitudes. This leads, through a two-step 
iterative procedure, to a two-loop representation 
of the K TTTTTT amplitudes. 

In the following we illustrate the method on the 
simplest case of the Kl Stt" decay, mainly be- 
cause the analytical expressions are less involved, 
but also because, for the time being, we want to 
avoid addressing some further issues appearing in 
the treatment of the processes involving charged 
pions. The cusp in this decay has been observed 
by KTeV [7], as well as by NA48/2 [8]. 

2. Reconstruction theorem 

The approach we wish to implement proceeds 
in parallel with the construction of the two-loop 
representation for the tttt scattering amplitude 
achieved in [9110] . We shall use the fact that 
it can be extended to other processes and that 
isospin symmetry is not an essential ingredient 
jllj . Here we are interested in Kn tttt, related 
to K 3tt by crossing symmetry. The essential 
ingredients required in order to implement this 
construction are the following. First, we need a 
decomposition of the amplitude of the type 

A{s,t,u) = 167r(/o(s) + 3/i(s)cos0)+A>2 (2) 

with the following chiral behaviour, 

Re^f>2~0(/), ImA>2~0(/), (3) 
Re/,-0(p2), Im/,^0(p4), £ = 0,1. (4) 

Having this, one can reconstruct the amplitude 
A{s,t,u) of a process AB — > CD to 0{p^) from 
the knowledge of the imaginary parts of S and P 
partial waves of all the crossed amplitudes: 

A{s,t,u) = P{s,t;u) + $o(s) 
+ [s(i-u) + (m^-m|)(m^-TO|,)]*i(s) 

+ crossed channels + 0{p^), (5) 

where P(s,i;w) is a polynomial having the same 
s, t, u symmetries as the amplitude A{s, t; u) and 
of at most third order in the Mandelstam vari- 
ables. $0 and <I>i are the dispersive integrals of 



the partial waves of the s-channel amplitude 

*.(.) = m.f (6) 

J thresh. ^-^ 

48.3 r d. 

Jthrcsh. X^{x — s)Xjlg{x)X(jjj{x) 

and similar for the t- and w-crossed channel 
[XAsix) = {s- {tua + rnBY){s - {tua - mB)^)]- 
The imaginary parts that enter the above ex- 
pressions are obtained from the unitarity relation, 

Im/r^(.).^l^ 

k 

x/r'W(//^'(5))*e(s-thr,), (7) 

projected on the corresponding partial waves. 
The sum goes over all the possible intermediate 
states fc (5 is a symmetry factor, S' = 2 for undis- 
tinguishable states and S = \ otherwise). In the 
low-energy region, and up to two-loops, these are 
restricted to pairs of light pseudoscalar mesons. 
To the extent that we are interested in the de- 
cay region only, we may further restrict them to 
intermediate tttt states. The contributions from 
other intermediate states, like e.g. Kit, can be 
expanded in powers of the Mandelstam variables 
and absorbed into the polynomial P(s,t,u). 

This unitarity relation and the reconstruction 
theorem can be used iteratively, i. e. starting from 
the LO amplitudes, we obtain the NLO results. S 
and P partial wave projections thereof then allow 
to obtain the NNLO expressions. Details of the 
first step will be presented in the next section. 

According to the reconstruction theorem and 
due to the crossing symmetry, the two-loop rep- 
resentation of the Kl ^ 37r° amplitude looks like: 

Al,oo{s, t, u) = Pl;oo + 1>o ^""(s) + K'^'^it) 

+ $o^^"°H + 0(/) (8) 

with the polynomial 

Pl;oo = Cf{A^,M], + {Co^„[(s - s^)'] 

+ E^„[{s - 4)3]} +{s^t} + {s^ u}), (9) 

where the centre of Dalitz plot was defined as 
Sq = 1/3M|- -I- m§ and Cp corresponds to the 
standard normalization, Cp ~ ~f ^JsKid^ • 
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3. First iteration: one-loop expressions 

As already stated, we need the leading order tttt 
and Kit — > tttt scattering amplitudes. From the 
chiral perturbation theory we know that at 0{p^) 
they are represented by polynomials of at most 
first order in the Mandelstam variables. Their 
particular choice (connected also with the partic- 
ular choice of the polynomial of the reconstruc- 
tion theorem) is important since different choices 
can possibly lead to different convergence prop- 
erties of the chiral expansion and affect the sta- 
bility of the fit to the data. The standard choice 
(believed to be stable) is the expansion in sub- 
threshold parameters. For the tttt amplitude, this 
corresponds to 
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Another possibility would consist in choosing the 
scattering length and effective range parameter 
as independent coefficient. It is even possible to 
adjust the polynomial part P(s,t;u) of the tttt 
amplitude so that these coefficients retain their 
physical interpretation up to two loops, just as in 
the non-relativistic approach [5]. We shall study 
this option elsewhere [T5]. In the case of i^^Tr^ — > 
tt'^t:^ and K^tt^ tt^tt^ amplitudes we have 

•^L- + - 

where 

The 0{p^) chiral perturbation theory result is 
reproduced by special values of the parameters 
(aoo = 1, /3±o = 1, a±o = {2m\_ ~ ml) /ml; cf. 
[13] ) and similarly for the Kn part (see [T4]). 

Using these amplitudes in the first iteration, we 
obtain the one-loop result for <I>q of ([5]) as 

*o [s) ^ —AQ^MKaQomQ.h{s) 



= CF[Bl_{s-si) + A^^_Ml], (11) 
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polynomial + O(p^). 
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Figure 1. The partial decay rate Kl — > 37r° 
(in arbitrary units) as a function of the invariant 
mass of the tt^tt" pair squared for one particular 
choice of parameters. Within the small frame the 
cusp region is zoomed. 



The loop functions J are defined by 

JpQi^) = Thr dx^^^, (13) 
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which simplifies for both masses the same to 
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(14) 



Our NLO parameterization of the amplitude 
jl]) with $0 given by ^ can be now fitted to 
experimental data, and so we can get the 0(jP) 
values of all the seven constants appearing in the 
LO parameterization of ([TU)) and pT|) . It is also 
important to have in mind that our parameteriza- 
tion of the amplitude encompasses the chiral per- 
turbation theory result [M] as a particular case 
for special values of the parameters. 

As an illustration of the cusp effect in the 
Kl 37r° decay, we plot the resulting partial 
decay rate for one particular choice of the sub- 
threshold parameters in the Fig. 1. 

4. Second iteration: some remarks 

The second iteration leads to the two-loop ex- 
pression of the K — > 37r amplitude. However, a 
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few complications have to be dealt with. As we 
have seen above, in order to get the discontinu- 
ities of the fimction <I>q''"^(x), we need an appro- 
priate analytic continuation of the unitarity rela- 
tion (and therefore also an analytic continuation 
of the S and P partial wave projections of the 
amplitudes obtained by the first iteration) below 
the physical threshold. In the isospin limit the 
solution is well known. The general method de- 
veloped in [T3] is based on the careful definition 
of the i-integration contour C{t-\.,t-) in the for- 
mula for the partial wave projections, where the 
endpoints are 

t±{s) - I (3,s^ - s ± Xli\s)ao) + ie, (15) 
• 9t±{s) 

signer sign (16) 

The contour C(i-|-,t_) is defined in such a way 
that it avoids the intersection with the cuts at- 
tached to branch points corresponding to the nor- 
mal threshold of the amplitude AL:Oois,t,u) in 
the t— and u— channel. The trajectory of t±{s) 
and the basic types of the contour C(t+,t_) in 
the complex plane are depicted in [15| . In the 
case of the K^-k^ tt'^tt'^ and Kltt'^ tt+tt^ 
scattering, the generalization of this prescription 
beyond the isospin limit is straightforward. How- 
ever, this is not the case of e.g. the experimen- 
tally more interesting process K'^ir"' — s- tt'^tt*', 
where the naive application of the prescription 
[15j shows some problems. Namely, for the re- 
construction of the amplitude of this process be- 
yond NLO we need to compute the partial waves 
of K^-K^ TT^TT^ , where the trajectory of t_(s) 
crosses the t— channel cut of this amplitude in- 
stead of avoiding it, contrary to the isospin limit. 
In addition, the two-loop amplitude for K^n^ 
tt^tt'^ suffers from (complex) anomalous threshold 
stemming from the triangle Landau singularity. 
The consequence is that a careful analytic contin- 
uation of the normal dispersive integrals entering 
the representation of the amplitude by means of 
the reconstruction theorem has to be performed. 
These technical issues deserve a detailed further 
discussion which is however beyond the scope of 
this talk. 



5. Conclusion 

The cusp effect in K ^ 3tt offers an interesting 
possibility to extract quantitative information on 
TTTT scattering lengths from the experimental data. 
We have outlined a construction, in a fully rela- 
tivistic framework; of the corresponding two-loop 
amplitudes, based on general properties, unitar- 
ity, analyticity, crossing, and chiral counting for 
the partial waves. Our analysis provides two pa- 
rameterizations, one in terms of the subthreshold 
parameters and the other, as in existing analyses, 
directly in terms of scattering lengths. 
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